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Abstract 

We performed a self-consistent calculation using T-Matrix approximation for a quasi-bidimcnsional 
system. We calculated the one particle spectrum function A (k, uj) in the presence of strong d-wave attrac- 
tive interaction. The c-axis charge dynamics was studied by considering incoherent interlayer hopping and 
afr-plane charge dynamics was studied in the coherent limit. It is shown that the c-axis charge dynamics 
is mainly governed by the scattering from the in plane fluctuations. We also present results for c-axis and 
afr-plane resistivity and for thermopower coefficient. 
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1 Introduction 

In recent years one of the significant advances in the research of high-T c superconductors is the detailed ex- 
periments on the pseudogap behavior in underdoped high-T c superconductors. The pseudogap phenomena 
have been recognized as one of the most important issues. There are enormous studies from both experimen- 
tal and theoretical point of view. Influence of the pseudogap formation shows up in different experimental 
probes, such as NMR relaxation time, Knight shift, neutron scattering, tunneling, photoemission, specific 
heat, optical conductivity and d.c. resistivity 

The pseudogap phenomena means the suppression of the low frequency spectral weight without any long 
range order. For cuprate superconductors the low frequency spectral weight begins to be strongly suppressed 
below some characteristic temperature T* much higher than T c . Typically T* is several times larger that 
T c in underdoped cuprates and their doping dependence is qualitatively different. While T c deceases with 
underdoping, T* increases in contrast. This suggest that T* does not follow the real T c line but instead 
some kind of mean-field critical temperature Tmf- 

In reality the origin of pairing interaction should be considered to be the antiferromagnetic spin fluctu- 
ations. There are studies dealing with pairing correlations obtained by the spin fluctuations on the basis of 
the fluctuation-exchange (FLEX) approximation ||. To explain the pseudogap phenomena many theories 
have been proposed, such as the argument based on BCS-BE cross-over ||, spin-charge separation induced 
by stripe order RVB theory fla], strong superconducting phase fluctuations |], 0], magnetic scenario near 
the antiferromagnetic instability [||, ^, pi]] and formation of charge density waves (CDW) or spin density 
waves (SDW) [[0]]. A feature of the last scenario is that the pseudogap can persist in the superconducting 
state The effect of the magnetic field on the pseudogap in high-T c cuprates has been theoretically 

investigated by Yanase et.al. |jl3|| considering the Landau quantization for the superconducting fluctuations 
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as the main effect of the magnetic field. The pseudogap phenomena was studied in the context of non- Fermi 
Anderson model in 14] considering the electron-fluctuation interaction. Recent dimensional crossover study 
by Preosti et.al. [15| shows that the pseudogap effect is basically absent in three dimensions. 

In this paper we describe the pairing scenario based on the strong coupling superconductivity. Generally 
the strong coupling superconductivity indicates the existence of the incoherent Cooper pairs as the Nozieres 
and Schmitt-Rink formalism [16|. We think of the pseudogap as the phenomena brought about by the strong 
superconducting fluctuations. Maly et.al. [17] have introduced the idea of "resonances" which means the 
scattering mechanism by the presence of metastable preformed pairs. 

Through several years of extensive experimental work general consensus regarding the superconducting 
gap symmetry seems to be reach, that the superconducting gap has mainly ci-wave character with possibility 
of a small mixture of other angular momentum state [18, 19, 20]. A remarkable point of the pseudogap is 
that its structure in momentum space is the same as the superconducting d x 2_ y 2 symmetry with continuous 
evolution through T c . This implies that the pseudogap phenomena has close connection to the supercon- 
ducting fluctuations pi]] . On the other hand using Intrinsic Tunneling Spectroscopy, Krasnov et.al. p2[| 
found that the pseudogap is coexisting with the superconducting gap which imply that the two phenomena 
are different in nature. The precursor superconductivity scenario can not explain these experimental data 
and a possible explication might be the formation of CDW or SDW [11|. Another striking characteristic of 
cuprate superconductors is the anisotropy of their structure. The quantity which most evidently displays the 
anisotropy of a particular material is the optical conductivity a (co,T) and the corresponding d.c. resistivity 
p(T) ||. 

The opening of the pseudogap has drastic effect on the physical properties of the high-T c cuprates [GJ. 
It is found that associated with the pseudogap the in-plane resistivity deviates from the T-linear behavior 
[24 1 and the T coefficient of the c-axis resistivity changes sign, showing semiconducting behavior at low 
temperature [25]. 

This paper is constructed as follow. In section 2 we give the model Hamiltonian and explain the the- 
oretical framework. In section 3 we present the results obtained for the spectral function A (k, to) and the 
results for the transport properties studied in this paper. In section 4 we summarize the obtained results 
and give conclusions. 



2 Theoretical Framework 
2.1 Model Hamiltonian 

In this section we describe the theoretical framework of this paper. We consider a microscopic model which 
incorporates both strong electron-fluctuation interaction in the CuO planes and a weak interlayer coupling. 
Within each layer we consider the following two-dimensional model Hamiltonian which has d x 2 _ y 2-w&ve 
superconducting ground state: 

Hi = £kC^ CT Q k(7 + ^kk' c ^T Qk +qt c «k'l Qk '- c ii W 

k,cr qi k,k' 

where Qko- is the annihilation operator for an electron with momentum k and spin a in the layer I. The 
electron dispersion relation is expressed as: 

£k = — 2t (cos k x + cos k y ) — 4t' cos k x cos k y — fi (2) 

where t and t' are the nearest and next-nearest neighbors hopping amplitudes and \x is the chemical potential. 
From now on we will consider t equal to unity and t' = — 0.5i. The pair interaction is written as: 

^kk' = Vhh' (3) 
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where 

/k = cos k x — cos k y (4) 

is the d x 2_ y 2-w&ve form factor. In the previous relation V is negative. The hopping between layers is 
included in: 

H = J2 H l ~ *c J2 ( C tia c l+U<T + Ci~+lia c Ucr) ( 5 ) 
I I 

where t c is the hopping amplitude between layers and cn a is the annihilation operator for an electron within 
the planar site i, with spin a in the layer I. The important character of high-T c cuprates is the momentum 



dependence of the interlayer hopping matrix element t c (k). The band calculations [gg, 27] of these materials 
have shown that the c-axis tunneling matrix element is larger at (jr, 0) and symmetry related points. The 
functional form is: 

t c (k) = (cos k x — cos k y ) 2 (6) 

On the other side the coupling between layers has little effect on the pseudogap p8| . The out of plane 
resistance is determined mainly by the properties of individuals CuO layers. The large magnitude of the 
resistivity anisotropy reflects that the c-axis mean free path is shorter than the interlayer distance, and the 
carriers are tightly confined to the CuO planes and also is the evidence of the incoherent charge dynamics 
in the c-axis direction. 

The interlayer part of the Hamiltonian (|5|) will be considered for the calculation of c-axis charge dynamics 
and p c (T). 

2.2 Self-Consistent T-Matrix Approximation 

In this paper we focus on the calculation of the Green function directly on the real- frequency axis [29] in 
order to avoid the difficulties of controlling the accuracy of the calculations used in the numerical analytical 



continuations from the imaginary to real axis by Pade algorithm [ 30 1 . 



We carry out a self-consistent calculation for the spectral function: 



A(k,u) = — ImG(k,u) (7) 
where the retarded Green function G (k, u) is given by: 

G(k,w) = 1 (8) 

UJ — £k — 2j (k, LO) 

where S (k, u) is the retarded self-energy. The imaginary part of the self-energy can be expressed as: 

JmS (k, lo)= flJ2 [ duJ ' if + n ( u + A ( k '> w ImT ( k + k', w + (J) (9) 
k' •* 

where / (u>) and n (u) represent Fermi-Dirac and Bose-Einstein distributions. The real part of the retarded 
self-energy can be calculated directly or using Kramers-Kronig relation: 

Reljik., uj) =p.v. / ; — (10) 

J 7T UJ — Uj' 

where p.v . represents the principal value of the integral. The Hartree term is neglected in the self-energy. 
The T-Matrix is given by the following relation: 

V 2 U(k,uj) 

T(k '" ) = TTW(k^) 
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where II (k, uj) can be evaluated directly from the spectral function as: 

ImU (q, uj) = vr J du'fk tanh ^ A ( k , u') A (q - k, 



UJ — UJ 



The chemical potential is determined by fixing the density n through the relation: 



n 
2 



J2 J duA(k,uj)f(co) 



(12) 



(13) 



We can calculate the hole doping 5 using the relation 5 = 1 — n. Beside relation ([13]) we must take care of 
the sum-rule satisfied by the spectral function: 



1 = J2j d ^A(k,uj) 



(14) 



which must be evaluated at each step in the self-consistent calculation. In the present calculation we divided 
first Brilouinne zone into 32 x 32 or 64 x 64 lattice and the frequency integration was replaced by a discrete 
sum of 1024 or 512 points. We did not found differences in the behavior of spectral function using these two 
discretizations. 

Performing the self-consistent calculation given by Eqs. (fz|)-(|l4|) we obtain the spectral function A (k, uj). 
We also calculated the density of states: 



N(uj)=Y,A(k,uj) 
k 



(15) 



Having obtained the spectral function we can calculate different physical properties of the cuprates. The 
real part of the aft-plane optical conductivity is given by: 



f(uj)-f(uj + u) ^ 

Cab [y) = 00 / duj 2^ 



\dk x J 



dk„ 



A(k,uj) A(k,uj + v) 



(16) 



where o"o = e 2 7r/2/i. 

The in-plane conductivity a a b (uj) given by eq. ( |l6|) neglects the vertex correction. The in-plane conduc- 
tivity is coherent in character and shows a Drude peak at low frequencies. 

For the calculation of c-axis conductivity a c (uj) we consider the incoherent case. Incoherent conductivity 
corresponds to diffusive c-axis transmission and amounts in taking the averages of the spectral function N (uj) 
over all momenta. 

/ (uj) - f (uj + v) 



a c (u) = a duj 



J2tl(k)N(uj)N(uj + 



(17) 



The same relation appear in the calculation of optical conductivity in disordered systems [pi] ]. We calculate in 
the limit uj — > the in-plane resistivity p a ^ = lim a; ^o(l/ "afe and c-axis resistivity p c = lim £t) _ > o(l / o c (uj)). 
In order to evaluate the thermopower we have to calculate p2^]: 



L 



'.J 



duj 



duj 



and the thermopower can be calculated using: 



S 



lei T Ln 



uJ 



3-1 



(18) 



(19) 
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3 Results 



3.1 Spectral Function and Phase Diagram 

In this section we present the obtained results for spectral function, using the self-consistent T-Matrix 
approximation described in previous section. In order to make the results more transparent we present 
the lowest order calculation results and the self-consistent calculations. We also present the phase diagram 
obtained for this model. As we can see both in Fig. 1 and Fig. 2 in the lowest order calculation we obtain a 
pseudogap along the direction (0, tt) — (tt, tt) and no pseudogap is obtained along the (0, 0) — (tt, tt) direction. 
This is in agreement with ARPES measurement s p3[ . Janko et.al. claims that this anomaly is caused by 
nearly stable cooper pairs in the excited states |34j]. The same behavior is observed in the self-consistent 
calculation as we can see in Fig. 3 and Fig. 4. Qualitatively we found small differences in behavior between 
lowest and self-consistent calculation for the spectral function. 

We also calculated the critical temperature T c of the system using Thouless criterion [lq| . When 1 + 
FIT (k, uj) = the superconductivity occurs. Thouless criterion is equivalent to BCS theory in the weak 
coupling limit. In our approximation T(k, uj) can be regarded as the propagator of the Cooper pairs. 
Thouless criterion correspond to the situation in which T(k, uj) has its pole at k = and u = 0. We are 
interested in the behavior of the system in the normal state, near the superconducting critical point, where 
the superconducting fluctuations are inhaced. 

In Fig. 5 we present the phase diagram of the model where T c and Tmf are the critical temperature 
calculate in the self-consistent and lowest-order calculated using Thouless criterion. We also calculated the 
critical temperature in the first order of perturbation theory and found that in the limit of strong coupling 
the self-consistent calculation of Tc reach the limit of first order calculation of Tc- Tmf can be viewed 
as the BCS critical temperature. The phase diagram corresponds to the underdoped region with the hole 
concentration 5 = 0.1. Roughly speaking the pseudogap phenomena occurs in the region between T c and 
Tmf- This region is extremely wide in the strong coupling case. As we can see from the shape of spectral 
function A (k, uj) the main effect of the electron-fluctuation interaction is the change of shape of A (k, uj) 
which is seen near k = (0, tt). This suggests that the electron-fluctuation interaction shifts the Fermi surface 
to the place far from k = (0,7r). Similar with other reported self-consistent treatments of the T-matrix 



J35| , 36] the pseudogap is not seen clearly in the integrated density of states but only in the k-dependent 
spectral function at extreme parameters ranges. In Fig. 6a we present the density of states for different 
couplings at the temperature T = 0.35t and in Fig. 6b the temperature dependence of the density of states 
for fixed coupling strength V = — 3t. As we can see from Fig. 6a, increasing the coupling strength, the 
density of states near u = slowly decreases but the central peak still remain. The central peak is due to 
the integration over the direction (0, 0) — (tt, tt) in the spectral function where no pseudogap appears. For 
the same temperature the central peak remain at the same frequency uj ~ 0. Decreasing the temperature 
the central peak start to shifts to lower frequencies and a hump start to appear in the density of states. The 
two humps appearing in Fig. 6a are due to the pseudogap formation in the spectral function near the (0, tt) 
and symmetry related points and their intensity increase with lowering the temperature. 

3.2 Optical Conductivity and Resistivity 

We have calculated the in-plane and c-axis optical conductivity using eqs. (|i~6|) and ([l7|). We also have 
calculated the resistivity as function of temperature for both in-plane and c-axis direction. Our approx- 
imation assumes the independent electron propagation in each layer and is justified for small interlayer 
hopping. The interlayer hopping hopping matrix element t c is weak resulting in an incoherent c-axis trans- 
port. Many approaches have been suggested [37, 38, 39] to describe the c-axis transport properties. The 



optical conductivity was calculated combining well controlled analytical and numerical methods in the case 
of one-dimensional Mott-Hubbard insulator at zero temperature in . The c-axis optical conductivity and 
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d.c. resistivity was calculated within t — J model in [41]. 



In Fig. 7 we present the obtained results for the in-plane conductivity as function of frequency for two 
temperatures. The in-plane spectrum of a a b (uj) is dominated by Drude peak at uj = 0. The intensity of the 
peak grows by increasing the temperature. In Fig. 8 we present the c-axis optical conductivity for three 
temperatures at the doping level 5 = 0.1. The electronic contribution to a c (uj) in the normal state of the 
underdoped cuprates is nearly uj independent and does nor form a Drude peak at u = 0, indicating incoherent 



nature of the charge dynamics along the c-axis [42]. It should be emphasize that the pseudogap appears 
only in a c (uj). Actually the development of the pseudogap in a c (uj) is correlated with a rapid narrowing of 
Drude width in a c (uj). We observe an anomalous broad peak in a c (uj) at finite frequencies which grows in 
intensity as the temperature is lowered with no discontinuity at the superconducting transition similar to the 
experimental results [43]. In Fig. 9 we present the temperature dependence of the in-plane resistivity p^. 
We found a linear behavior for the in-plane resistivity at high temperature and a deviation from linearity 



below a characteristic temperature T*. T* may be defined by different methods [44]. In this paper we 
consider the temperature T* the characteristic temperature where the resistivity present a downturn, with 
a linear behavior above it. The arrows in Fig. 9 indicate T* . Increasing the doping concentration 8, T* 
decrease. 

In Fig. 10 we present the c-axis resistivity. We found in this case the same linear behavior of resistivity as 
function of temperatures, at high temperatures and a upturn at lower temperatures (dp c /dT < 0) signaling 
a semiconductor -like behavior. The semiconducting T dependence of the c-axis resistivity is associated with 
the deepening of the pseudogap. Therefore the contribution to the conductivity from the quasiparticles near 
the (tt, 0) is suppressed by the pseudogap. The in-plane conductivity is determined by the contributions from 
"cold-spots" and the pseudogap is not important for the in-plane conductivity, but has more effect on the c- 
axis conductivity. The c-axis conductivity is determined from the contributions from the "hot-spots" where 
the pseudogap is large. Therefore the c-axis conductivity is remarkably suppressed by the pseudogap. This 



is in agreement with the experimental results for cuprates [45|. Our results are qualitatively in agreement 



with the experimental results obtained for (R\_ x Ca x ) Ba2Cu^OQ + s where with Ca doping the number of 



holes in CuOi planes increases p6|, 47], where R stands for Y and Lu. p a b is linear in T above T* and 
decreases more steeply below T* |24[| . 

3.3 Thermopower 

Thermopower has become a simple and ideal probe for the studies of high Tq superconductors. The ther- 
mopower is one of the several physical quantities which distinctly reveals the unusual normal state properties. 



The results obtained for the thermopower are calculated using eq. (|19[). The thermopower of cuprates ||48[] 
has the propriety that is not a monotonic function of temperature. Several authors have tried to explain 
the temperature behavior of thermopower in term of phonon drag, metallic diffusion, semiconducting model 
and others p^], but they cannot explain the linear behavior at high temperature [^] without considering 
the cylindrical Fermi surface pi]]. 

In Fig. 11 we present the temperature dependence of the thermopower. It can be seen that the normal 
state thermopower decrease with increasing hole doping 5 to a lower positive value, at a fixed temperature. 
The temperature dependence of thermopower is linear with a negative slope at high temperature and present 



a broad peak, shifted toward the low-temperature as the doping increase [52]. The hump of thermopower 
in the normal state was reported as a kind of universal behavior of HTSC, corelated with the opening of 
the pseudogap in the normal state of a superconductor [m ]. 

In Fig. 12 we present the doping dependence of the thermopower at the interaction strength V = — 5.5i 
and at the temperature T = 2t. At a fixed temperature much larger than the critical temperature, S 
continuously decreases from underdoped to overdoped regions, changing its sign near optimal doping. The 
theoretical work presented in j54| shows that the doping dependence of the thermopower can be explained 
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by the common band dispersion relation of high-T c cuprates. Therefore the thermopower can be a criterium 
for measuring the hole doping in CuO layers [55|. 



4 Conclusions 

In this paper we have shown that the d-wave pairing scenario, based on the strong coupling superconductivity 
well explains some of the physical properties of the layered cuprates above the critical temperature T c . We 
have considered that the strong coupling superconductivity takes place in underdoped cuprates and the 
pseudogap phenomena are its precursor. 

We considered the T-Matrix approximation which is the dominant scattering vertex in the vicinity of 
the superconducting critical point. In this analysis we didn't considered the vertex corrections. Anyway 
our theory appropriately explains the pseudogap phenomena in high-T c cuprates. We did not take care 
into account the antiferromagnetic spin fluctuations which are important for explaining other properties of 
cuprates such as the NMR spin lattice relaxation \/T\. 

Based on this scenario we studied charge dynamics, in-plane and c-axis resistivity and thermopower and 
compared the results with the experimental results. As can be seen from previous section our theory fit well 
the experimental results from YBCO class of layered superconductors. For the calculation of in-plane charge 
dynamics we obtained a spectrum dominated by a Drude term, while for c-axis dynamics we found an almost 
frequency independent behavior. For the c-axis resistivity we found, due to the opening of the pseudogap 
in the density of states, a semiconductor-like behavior. For the temperature dependence of thermopower 
we found a linearly dependence at high temperatures with a negative slope. The negative slope depends 
weakly on the doping level. 
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Figure Captions 



Figure 1 

Spectral function in the lowest order approximation at T = 0.6t in the (0, ir) — (ir, ir) direction at the 
doping concentration 5 = 0.1 and with the coupling strength V = — 5t 
Figure 2 

Spectral function in the lowest order approximation at T = 0.6t in the (0, 0) — (ir, ir) direction at the 
doping concentration 5 = 0.1 and with the coupling strength V = — 5t 
Figure 3 

Spectral function in the self consistent approximation at T = 0.6t in the (0, ir) — (ir, ir) direction at the 
doping concentration 5 = 0.1 and with the coupling strength V = — 5t 
Figure 4 

Spectral function in the self consistent approximation at T = O.Qt in the (0, 0) — (ir, ir) direction at the 
doping concentration 5 = 0.1 and with the coupling strength V = — 5i 
Figure 5 

Phase diagram of the model at hole concentration 8 = 0.1 
Figure 6 a 

Density of states for temperature T = 0.35i and for different coupling strength. The solid line correspond 
to V = —t, the dotted line to V = —3t and the dashed line to V = — 6i. The hole concentration is 5 = 0.1. 
Figure 6 b 

Density of states for different temperatures and for coupling strength V = —3t. The solid line correspond 
to T = 3.5t and the dashed line to T = 0.35t at hole concentration is 5 = 0.1. 
Figure 7 

In-plane optical conductivity at T = OAt (solid line) and T = 3At (dashed line) at the doping concen- 
tration 5 = 0.1 and coupling strength V = —5.5t 
Figure 8 

c-axis optical conductivity at different temperatures. T = 0.6t (solid line) T = l.Ot (dash line) T = l.6t 
(dot line) in the limit of small frequencies at the coupling strength V = — 5.5i and doping level 5 = 0.1 
Figure 9 

Temperature dependence of the in-plane resistivity at different densities, 5 = 0.1 (solid line) (Tq = 0.21t), 
5 = 0.15 ( dashed line) (Tc = 0.3£) and 5 = 0.2 (dashed-dotted line) (Tc = 0.36i), at the same coupling 
strength V = —5.5t. The arrows indicate the T* . 

Figure 10 

Temperature dependence of the c-axis resistivity at different densities, 5 = 0.1 (solid line) (Tc = 0.21i), 
5 = 0.15 ( dashed line) (Tc = 0M) and 5 = 0.2 (dotted line) (Tc = 0.36t), at the same coupling strength 
V = -5.5*. 

Figure 11 

Temperature dependence of thermopower at different densities, 5 = 0.1 (solid line) (Tc = 0.21t), 5 = 0.15 
( dashed line) (Tc = 0.3t) and 5 = 0.2 (dotted line) (Tc = 0.36*), at the same coupling strength V = — 5.5t. 
Figure 12 

Termopower as function of the doping concentration 5 at the coupling strength V = —5.5t and T = 2t 
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